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Entanglement from density measurements: analytical density-functional for the 
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We derive an analytical density functional for the single-site entanglement of the one-dimensional 
homogeneous Hubbard model, by means of an approximation to the linear entropy. We show that 
this very simple density functional reproduces quantitatively the exact results. We then use this 
functional as input for a local density approximation to the single-site entanglement of inhomoge- 
neous systems. We illustrate the power of this approach in a harmonically confined system, which 
could simulate recent experiments with ultracold atoms in optical lattices as well as in a superlattice 
and in an impurity system. The impressive quantitative agreement with numerical calculations - 
which includes reproducing subtle signatures of the particle density stages - shows that our density- 
functional can provide entanglement calculations for actual experiments via density measurements. 
Next we use our functional to calculate the entanglement in disordered systems. We find that, at 
contrast with the expectation that disorder destroys the entanglement, there exist regimes for which 
the entanglement remains almost unaffected by the presence of disordered impurities. 
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I. INTRODUCTION 

Entanglement, one of the most intriguing features of 
quantum mechanics, is an essential ingredient in quan- 
tum information theory (l|, and it is also receiving in- 
creasing attention in the study of fundamental phenom- 
ena such as quantum phase transitions (QPT) |2J . Entan- 
glement has been investigated in many different physical 
systems, such as photons in optical cavities Q, ultracold 
atoms interacting with light [J, |5[ , and in solid state sys- 
tems [2[ , where it has been associated to particle-particle 
correlations Q. 

Recently, special attention has been devoted to solids, 
since they are considered one of the most promising sys- 
tems for the development of quantum information devices 
Q. In particular, the Hubbard model j8[, which cap- 
tures most important physics of interacting many-body 
systems |9HTa|. has been intensively used to understand 
entanglement in solids [13HL9J . In this context, creat- 
ing, distributing, quantifying, and controlling entangle- 
ment represent great challenges. Much progress has been 
achieved and some well defined entanglement measures, 
such as entropies and concurrence, have been proposed 
[2J]. The von Neumann entropy has been used for mea- 
suring the bipartite entanglement of pure states [2j . For 
systems with very many degrees of freedom, the linear 
entropy - which gives an indication of the number and 
spread of terms in the Schmidt decomposition of the state 
- has been proposed as a viable alternative |20l - [23 |. 



For the homogeneous Hubbard model, the single-site 
entanglement |23[ and its properties have been investi- 
gated 1H and associated to QPT's QJ, QUI]. A local- 
density approximation (LDA) for the entanglement en- 
tropy of any generic inhomogeneous system has been pro- 
posed [lj| and applied to the Hubbard model for different 
types of inhomogeneities. Also block-block entanglement 
has been investigated in connection to QPT [25[ and from 
the view point of its universal and non-universal contri- 
butions [261 ] . 

Another extremely important - but less addressed - 
question concerning entanglement quantifiers, is whether 
and how theoretical predictions can be connected to 
quantities that are measured experimentally. It has been 
shown that magnetic susceptibility can be an entangle- 
ment witness |27H29| , which is used to point out the pres- 
ence of entanglement. Also momentum-momentum cor- 
relation [30, [3jJ, energies [2ll |22|. and probability den- 
sity at highly symmetric points [221 ] have been proposed 
as entanglement indicators, which reproduce the general 
trend of entanglement. However, a quantity able to pre- 
dict entanglement quantitatively and, in addition, with 
the potential of being experimentally measurable, is still 
missing. Here we show that a very promising candidate 
for this is the particle density. 

As a direct consequence of the Hohenberg-Kohn theo- 
rem [32|, the pillar of Density Functional Theory (DFT), 
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ground-state entanglement must be a functional of the 
density. For the one-dimensional homogeneous Hubbard 
model this functional may be numerically accessed [33| . 
but an explicit expression cannot be obtained due to the 
implicit dependence of the ground-state energy on the 
density. Different analytical and numerical techniques 
have been proposed for this calculation [13l - [l5l ]. Nev- 



ertheless, obtaining numerically the ground-state energy 
has been so far the main limiting factor for a complete 
analytical understanding of the entanglement in the ho- 
mogeneous Hubbard model. For inhomogcncous systems 
the scenario is even more complicated, despite the appli- 
cation of powerful DFT methods [16[ . 

From the experimental point of view, it has been re- 
cently shown that the measured density in systems of 
ultracold atoms in optical lattices with single-site reso- 
lution [33, (38| can be very accurately reproduced by the 
single-site density of the Hubbard model. Several the- 
oretically predicted phenomena have already been ver- 
ified in such experiments, e.g., the supcrfluid to Mott- 
insulator transition [39| and BEC-BCS crossovers j40| ; 
via DFT, the density may now become a powerful tool 
for measuring entanglement in experiments. 

Motivated by this, in the present work we provide an 
analytical expression for the density-functional for the 
entanglement of homogeneous strongly correlated sys- 
tems. Our approach is based on an approximation of the 
linear entropy for the one-dimensional Hubbard model 
valid, in principle, for a specific range of interactions and 
densities. The functional is then extended to a vast range 
of parameters by means of well established transforma- 
tions. We test the performance of our density functional 
in homogeneous and, via an LDA, in three different in- 
homogcncous systems. We find that the entanglement 
in many-body systems can be quantitatively obtained di- 
rectly from our very simple density functional. Finally we 
use our method for investigating the entanglement in dis- 
ordered systems. To achieve a good statistics, this anal- 
ysis requires many realizations and, therefore, could not 
be explored so far, while our functional greatly reduces 
its computational cost. We find that in general disorder 
destroys entanglement. Surprisingly though, there exist 
regimes for which the degree of entanglement remains 
almost unaffected. 



II. ANALYTIC DENSITY FUNCTIONAL FOR 
THE HOMOGENEOUS SYSTEM 

We start from the homogeneous Hubbard model, 
H = - X! (cJ CT Ci+i,ff + H.c.J + u ^ %"4 + V ^2 "iff) 

ia i ia 

(i) 

in which all sites are equivalent. Here u, the on-site inter- 
action, and V , the external potential, arc rescalcd by the 
hopping parameter between neighboring sites; c ia and Ci a 
are fermionic creation and annihilation operators, and 
hicr = &i a Ci a is the number operator at site i with spin 

o-=t,4» 

The single-site entanglement is quantified using the lin- 
ear entropy 

U = 1 - T^), (2) 

where the reduced density matrix of site i, pi = Trsp, 




FIG. 1: (Color online) Single-site entanglement as a func- 
tion of density for the homogeneous Hubbard model in the 
thermodynamic limit. For repulsive interactions (u = 2,4,8) 
the exact results (solid lines a, b, and c) are compared to 
Eq.® (dashed lines a, b, and c). For attractive interac- 
tions (u — —4) the exact results (solid line d) are compared 
to Eq. (|ll|l (dashed line d). Inset: zoom of main panel for 
n « a(u) + 1/2. The exact results were obtained from Eq.([3]|. 



is the trace of the total density matrix p with respect 
to the remaining sites. For the homogeneous Hubbard 
model pi = diag{w^,w\.,W2,wo}, with w a = (h a ) — w 2 
the single spin occupation probability, w% = (n-j-nj,) the 
double occupation, and wp = 1 — w^ — w^ — W2 the zero 
occupation probabilities |23| . 

This allows us to express the single-site entanglement 
of the homogeneous Hubbard model as an explicit func- 
tion of the occupation probabilities 



L hom = 1 - w 2 - u>l - w\ - w 2 . 



(3) 



For non-magnetic systems the probabilities are given by 



io 2 

'"0 



deo(n, u) 

du 
n 

= W ^2 
= 1 - n + w 2 , 



w 2 



(4) 

(5) 
(6) 



where eo(n,u) is the per-site ground-state energy. 

In order to express the entanglement as an explicit 
function of the density, we will approximate Eq. (QJ . Ini- 
tially, we focus on repulsive interaction and n < 1, but 
as shown later, the results will be extended to n > 1 
and to u < by means of particle-hole transformations. 
An exact analytical expression for the energy, and con- 
sequently for u>2, is unknown in this regime, except at 
n = 1, for which the exact Lieb-Wu solution [33j gives 



ui2{n = 1, u > 0) = a(u) 



J (x)J 1 (x)e ux / 2 



(1 



C /2)2 



dx, 



(7) 



where Jk{x) are Bessel functions of order k. 

For low-density systems the double occupancy is ex- 
tremely reduced due to the repulsive interactions |4l|; 
for higher densities, W2 increases and can not be dis- 
carded. We then consider W2 ~ when L hom (n,u > 
0)\ W2 =o > L hom (n,u > 0)| W2=a ( u ), and assume Eq.© as 
an approximation for u>2 otherwise [421 ] . 

Then the linear entropy can be written as the explicit 
density functional 

L hom {n,u>0) w 2n-^-+ (8) 

[{An - 2)a(u) - 4a(u) 2 ] Q(n - a(u) - 1/2), 

where Q(x) is a step function with 9(x) = for x < 
and <d(x) = 1 for x > 0. We extend our functional to 
the regime of n > 1 using the particle-hole symmetry, 
by replacing n with 2 — n in Eq.®. Such extension is 
necessary to describe a larger variety of systems, and also 
in view of extending the results via DFT methods to more 
realistic inhomogeneous systems. 

For attractive interactions the double occupancy can 
neither be discarded nor approximated by the Lieb-Wu 
result. The extension of our analysis to u < can be 
though obtained by means of the mapping between the 
energies of repulsive and attractive systems [44| , 

UTl / 7TTI \ 

e (n, u < 0) w — + e (n = 1, u > 0) sin ( — J , (9) 
which is valid for any density. From @ and © we obtain 



w 2 {n,u < 0) « - -a(|u| 



7TO\ 

TJ' 



(10) 



and the density-functional for the entanglement becomes 

,2 



- /torn 



(n, u < 0) « n 



_+2a(M)Bin(— J 



-4a 2 (|u|)sin 2 



2 /7m\ 



2 / 



(11) 



In Figure [T]wc compare our density- functionals for the 
entanglement, Eq.© and Eq. ([i"Tj) . to the exact entangle- 
ment Eq.([3|). The agreement is very good, considering 
the simplicity of our density-functional expressions. For 
u < the exact results are very accurately reproduced be- 
cause the double occupation is fully taken into account 
in Eq. fTTj) . For u > the results reveal that the rele- 
vant dependency on n for the entanglement is taken into 
account by the other occupation probabilities, i.e., the 
main contribution of the double occupancy to the entan- 
glement is almost independent on the density, therefore 
the simplified form of u^ adopted in Eq. © is enough for 
predicting the entanglement with good precision. In par- 
ticular Figure Q] shows that there are ranges of densities, 
< n < 0.5 and 1.5 < n < 2, for which many-body 
interactions do not play a relevant role: in a homoge- 
neous system, up to n w 0.5 (and for n > 1.5 when 
considering particle- hole transformation), the double oc- 
cupation is in fact strongly reduced as here W2 ^ ?i 2 /4 
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FIG. 2: Single-site entanglement as a function of the inhomo- 
geneity: a) k, the curvature of the harmonic confinement, 



Vi 



k(i - io) 2 , with j = N/2 + 1, N = 100, u 



and n — 0.7; the solid line is the derivative dL LD /dk from 
Eq. (|12[) with Eq.((8| as input and the arrows indicate stages 
of the density profile; b) AV, the modulation amplitude of a 
superlattice with periodicity of 20 lattice sites, for N = 200, 
100 particles and fixed interaction u — 3; c) Vi, the inten- 
sity of a single- impurity localized at the site i = N/2 + 1, for 
iV = 100, 80 particles and u — 6. 



(u>2 ^ (2 — n) 2 /4). In turn in the Hubbard model, where 
only on-site interactions are considered, this implies a 
lesser role of the particle-particle interactions [3J] . 

However there is a small range of combined parame- 
ters - namely n sa a(u) + 1/2 for < u < 2 (see inset 
of Fig. []} - for which our approximation to W2 is not 
as good: here also because on-site repulsion is small, u> 2 
may become significant. The maximum error in this spe- 
cific region decays exponentially from ~ 12% at u = 0.5 
to 4.6% at u = 2. Considering that typical interactions in 
solids correspond to u ~ 6 [35| and in experiments with 
optical lattices the interaction can be easily tuned in a 
wide range (up to u ~ 200) [36| , for most systems of prac- 
tical interest our approach can provide the entanglement 
with a precision better than ~ 3% for any density. 

The very good agreement of our density function- 
als with exact calculations shows that the entangle- 
ment of actual experiments described by the Hubbard 
model [33, |38[ could be obtained accurately and directly 
from density measurements. Also our results imply much 
simpler calculations which opens possibilities for the en- 
tanglement investigation in more complex systems, such 
as the disordered media we will explore in the last sec- 
tion. 



III. DENSITY FUNCTIONAL FOR 
INHOMOGENEOUS SYSTEMS 

In order to investigate the performance of our density 
functional in inhomogeneous systems, we present in what 
follows numerical calculations for the Hubbard model 
with an inhomogeneous external potential £^ a Vifii tiy . 
We adapt to this purpose the LDA formalism developed 
in Rcf . [16| . The average single-site linear entropy of in- 
homogeneous systems is then given by 



,v 



L mh (n,u)^L LDA (n,u) = - 



N ^' 



- horn 



(n,u)\ n=ni 



(12) 
In Eq. (fi"2"f the inputs are our density- functionals, Eqs.® 
or (TTTJ) , calculated at the local i-site density rii, and av- 
eraged over the N sites [43 |. 

We start with the inhomogeneous harmonic potential, 
which may model, for example, the trap in systems of 
ultracold atoms in optical lattices [37] or quantum dots 
[la . USh |45| . This inhomogeneous system contains very 
interesting physics: it undergoes subtle transitions in the 
density profile as the potential increases [46[ and the 
single-site entanglement is able to point out each and 
every one of these phases [Tg]. We have chosen n = 0.7, 
which represents a typical density reached in the experi- 
ments [471 ] . and u = 8, as ex per imentally the interaction 
can be tuned up to u ~ 200 [3j6J . In Figure [5]-a wc com- 
pare our density-functional for the entanglement to ac- 
curate numerical results. The agreement is indeed very 
good and all the stages of the particle density [16( are 
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FIG. 3: Upper panel: average single-site entanglement in dis- 
ordered chains as a function of the impurities' intensity V, 
for increasing concentrations C of randomly distributed im- 
purities, for 100 sites, 80 particles and u = 4. C is defined 
with respect to the total number of sites and each point is an 
average over 100 samples. Inset: entanglement as a function 
of V for C — 40% and different values of it. Lower panel: 
as for the upper panel, but for C — 20%, 80%, and 50%, to 
emphasize the C —¥ 100 — C symmetry for V — > — V of the 
entanglement curves. 



correctly pointed out by the derivative of our functional 
(see the entanglement derivative in FigJ^a). 

We have also checked the performance of our approach 
for different inhomogeneities and different values of u, 
such as a superlattice defined by a periodic modulation 
Ay in the external potential (u = 3), and the presence of 
a single impurity with intensity T^ = 5i (u = 6) as shown 
in Fig. [5]-b and Fig. [5]-c, respectively. The agreement is 
quantitatively similar to that found for the harmonic con- 
finement: for the superlattice the deviations are smaller 
than 2%, while for the impurity system we found devia- 
tions of about 2.2% for all impurity intensities (note the 



scale in the y-axis). 

Our results for three different physical systems con- 
firms that the entanglement can be calculated di- 
rectly from density measurements by using our density- 
functional even for inhomogeneous systems. As all the 
above tests have confirmed the reliability of our func- 
tional, we next use our density-functional for investigat- 
ing disordered systems. 



IV. ENTANGLEMENT IN DISORDERED 
SYSTEMS 

The investigation of any physical property in disor- 
dered systems usually requires an average over a large 
number of samples in order to improve statistics and 
avoid results peculiar to specific realizations. Thus, for 
example, for the generation of a single point in a system 
of 100 sites considering 100 samples, it would be neces- 
sary to perform 10 4 calculations of the desired average 
single-site property. This represents a serious limiting 
factor for estimating the entanglement via methods in- 
volving the computationally-expensive ground-state en- 
ergy calculation. In this sense our density-functional pre- 
scription represents a great advantage. 

Here we consider the average single-site entanglement 
over 100 realizations of a disordered system defined by a 
concentration C of randomly distributed impurities with 
intensity V. In Figure [3] we present the entanglement 
obtained with our density functional as a function of V. 
The entanglement is maximum for non disordered chains 
(V = 0) and for V < a plateau is observed for C < 70%. 
This plateau in the entanglement reflects the competition 
between the repulsive interaction u and the impurities' 
attractive potential V at the Mott-insulator transition 
induced by n% = 1 at the impurities' sites. A plateau is 
also observed in systems with a single impurity (see Fig. 
Hl-c), where the percentage drop of entanglement with the 
impurity strength is however much smaller. Interestingly, 
for C > 30%, we also find a plateau for V > (see Fig. 131 
upper panel). The physics behind both plateaus is simi- 
lar, but for repulsive V and u one needs a higher impurity 
concentration to reach the Mott-insulator transition, as 
in this case it is induced at the sites with no impurities. 
The reversed role between sites with and without im- 
purities manifests itself in the reflection symmetry with 
respect to V = for systems with C and C = 100 — C, 
as shown in Figj3]- lower panel. The width and height of 
the plateaus increase towards smaller impurity intensities 
with increasing (decreasing) impurity concentrations for 
V > (V < 0); the plateaus' width also increases with 



u as | V\ < u, as a result of the competition between the 
two potentials (see inset of Fig|3|). 

In disordered chains, similarly to what has been ob- 
served for other inhomogeneous systems, inhomogeneities 
are in general destructive to the entanglement, and the 
higher C and |V|, the more the entanglement is de- 
stroyed. However our calculations show that there are 
regimes for which the degree of entanglement is almost 
unaffected. The entanglement in this region is almost as 
higher as the homogeneous case, V = 0. This property 
could be useful for indicating the regime of parameters 
for which real-life devices may present a high degree of 
entanglement despite the presence of impurities. 



V. CONCLUSIONS 

In conclusion we have derived an analytical density 
functional for the entanglement in the Hubbard model 
represented by a very simple mathematical expression. 
We have shown that it reproduces the entanglement 
quantitatively in homogeneous and inhomogeneous sys- 
tems. We tested our functional on very different inho- 
mogeneous systems, i.e. an harmonic confinement, which 
simulates experiments of trapped atoms in optical lat- 
tices, as well as on a supcrlattice and a single-impurity 
system. 

We applied our density-functional prescription to the 
investigation of the entanglement in disordered systems: 
this analysis would be computationally very expensive - 
and as such difficult to achieve - with other methods. We 
found that, although in general the disorder destroys the 
entanglement, there exist a regime of parameters, both 
for repulsive and attractive impurity potential, for which 
the degree of entanglement remains almost unaffected. 
This result is promising as it suggests potential robust- 
ness for future device production. 

In addition of allowing simpler numerical calculations 
and paving the road to analytical analysis of many-body 
and quantum information properties, our findings are of 
direct interest for measuring the entanglement from ex- 
perimental setups. 
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